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An evolution equation for the expectation values of the Boltzmann factor between monomer 
valence bond states is derived. It contains the whole information on the thermodynamical and 
magnetic properties of the spin ^ quantum Heisenberg model. A new approximation scheme is 
proposed and studied in the thermodynamical limit. 
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0\ . I. INTRODUCTION 

In 1938 Hultheni determined eigenvalues and eigenstates of the spin \ quantum Heisenberg model on small clusters 
with V < 10 sites. At that time, computers were not yet at his disposal. He therefore exploited the symmetries of the 
<^ ' model-including the conservation of total spin S- in order to reduce the dimensions of the block Hamiltonians as -much 
ly-^ | as possible. With the advent of powerful computers exact diagonalizations on larger clusters up to V — 36 sitesa and 
- quantum Monte Carlo simulations on much larger systemsS could be performed. Many of these computations were 
devoted to thendetermination of the ground-state properties of the antiferromagnetic model: the energy and staggered 
CNI i magnetizationH, spin-spin correlatorsO etc. Since the ground-state is known to have spin OS it is sufficient to diagonalize 
the Hamiltonian in this sector. Conservation of total spin, however, is not used in most of these computations. An 
exception is Andersons resonating valence bond (RVB) stateQ. It seems that Hulthen's construction of spin states 
by means of valence bonds is not so suited for numerical calculations on large systems for the following reasons: 
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The set of valence bond states is overcomplete; there are linear dependences among them. 



• Valence bond states are not orthogonal to each other. 



There have been made some attempts to overcome these deficiencies!!. E.g. one can try to restrict the number 
of 'allowed' valence bond configurations in order to avoid over completeness. Such a procedure was successful for 
d = 1 u. For d > 1, however, it looks very unpleasant: The apparent symmetries of the system -translations, 
. rotations, reflections are hidden in complicated relations. 

In this paper I want to draw the attention to some astonishing properties of the projection operators on valence bond 
Q \ states: 

o 

• Completeness 

• Spectral like representations of the Hamiltonian and other operators. 

• Anticommutation relations with the Hamiltonian. 



The properties are discussed in section 2 and extended to projection operators on spin S states in section 3. In 
section 4 these properties are shown to describe the evolution of the expectation values of the Boltzmann factor 
between monomer valence bond states. In section 5 an approximation scheme for the evaluation of the evolution 
equation expansion is proposed: the monomer dimer cluster expansion. The thermodynamical limit of the lowest 
order approximation is investigated in section 6. 

II. VALENCE BOND PROJECTION OPERATORS (VBPO) 

We consider a lattice with an even number V = 2v of sites. Spin operators a{x) live at each site x. The operator 

(x,y) = ±(l-a(x)*(y)) (2.1) 
projects the two spin state at sites x and y onto its valence bond (i.e. spin ) part: 
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[x,y] = -j={x+(x)x-(y) - x-(x)x+ (</))• (2-2) 

Each partition of the lattice into pairs [xj,yj] of sites defines a valence bond configuration (VBC) K. We associate 
to each configuration K a state: 

\K) = [xi,yi\[x 2 ,y 2 }...[x v ,y v ], (2.3) 

with total spin 0. These valence bond states yield an overcomplete non orthogonal basis in the total spin sector. In 
the following we are concerned with the valence bond projection operators (VBPO): 

P(K) = \K)(K\ = (x 1 ,y 1 )(x 2 ,y 2 )...(x v ,y v ). (2.4) 

The VBPO have the following properties: 

tr P(K) = 1, (2.5) 



J2P(K) = a l , (2.6) 



K 

where lo is the projection operator onto the sector with total spin 0. The 'completness relation' (2.6) is proven in 
appendix A. If we take the trace of eq.(2.6) we get for the constant a : 

V" 

- " Mv-y 

where V\\ is the number of valence bond configurations and n (v) the number of spin states on a lattice with V = 2v 
sites. In appendix A it is shown that spin-spin correlators (a;, y) as they are defined in eq.(2.1) have a particular simple 
expansion in terms of the VBPO's: 

l (x, y)l = YAl V S([x, y] G K)P(K). (2.8) 
4a - 

According to eq.(2.8) all configurations K with a fixed valence bond [x,y] contribute with equal weight. The re- 
mainig configurations have zero weight. From eq.(2.8) we obtain the expansion of the nearest neighbour Heisenberg 
Hamiltonian projected onto the spin sector: 

H Q = 1 4(a;,y)l , (2.9) 

<x,y> 

ffo^^ftW), (2.10) 

K 

where 

N l (K)= J2 S([x,y]&K) (2.11) 

<x,y> 

is the number of dimers -these are the nearest neighbour valence bonds- on the configuration K. Eq.(2.10) looks like 
a spectral representation. But it is not, since the valence bond states \K) are not orthogonal. The Hamiltonian Ho 
obeys anticommutation relations with the VBPO's which are proven in appendix A as well: 

±[H ,P(K)] + = (dV + 3N 1 (K))P(K)+ (P(Q+(x,y)K)-P(Q_(x,y)K)). (2.12) 

<x,y> 

The sum on the right hand side extends over configurations which are obtained from the configuration K = 
[x,yi][x 2 ,y]--- by the following permutations: 

Q+{x,y)[x,yi][x2,y] = [x,y][x 2 ,yi], (2.13) 

Q_{x,y)[x,y{\[x 2 ,y] = [y, yi][x 2 ,x]. (2.14) 

Eqs.(2.5), (2.6) and (2.12) will provide us with an evolution equation of the Boltzmann factor between valence bond 
states. 



2 



III. PROJECTION OPERATORS IN SPIN S SECTORS 



For the construction of states with total spin S, we follow the original work of HulthenS: 

\Ks) = [z 1 ,z2,...z2s)s[xi,yi]...[x v ,y v ). (3.1) 

Here [z\, Z2, ■■■Z2s]s denotes a totally symmetrized state of spins at sites z%, ...Z2S- As before the [xj, yj] are valence 
bond (spin 0) states. Therefore we can characterize each of these states by a configuration Ks with v = V/2 — S 
valence bonds and 2S 'monomers'. Each site of the lattice is occupied cither by one valence bond or one monomer. 
The corresponding projection operators 

P(Ks) = \K s )(K s \, (3.2) 

have the following properties: 

tr P(K S ) = 1, (3.3) 

Y,P(K S ) = a s ls, (3.4) 

Ks 

where Is is the projection operator onto the sector with total spin S and: 

n s (y) 

fsiy) is the number of monomer- valence bond configurations Ks and ns{V) the number of spin S states on a lattice 
with V sites. The nearest neighbour Heisenberg Hamiltonian projected onto the spin S sector 

H s = ls ]T i(x,y)l s , (3.6) 

<x,y> 

can be expanded in terms of the projection operators P(Kg): 

Hs = ^-^y2N 1 (K s )P(K s ), (3.7) 

aS Ks 

where 

N t {K s )= S([x,y]eK s ) (3.8) 

<x,y> 

is the number of dimers on the configuration Ks- The Hamiltonian Hs obeys anticommutation relations with the 
projection operators P(Ks), which are proven in appendix A as well: 

hfls, P(K S )}+ = (dV + ZN^Ks) - 4N 2 (K S ))P(K S ) 

+ (P(Q+(x,y)K s ) - P(Q-(x,y)K s )). (3.9) 

<x,y> 

A^(-K's) is the number of monomer pairs on neighbouring sites < x,y >. The action of the permutation operators 
Q±<y{x,y) on the valence bond states is given in eqs.(2.13-14) and depicted in Fig. la). If monomers are involved, one 
has the following rules: 

Q+(x,z)[z,Z2, -} s [x,y} = [y,Z2,...] s [x,z] (3.10) 

Q-(x,z)[z,z 2 , ■■■} s [x,y) = [x,z 2 , ...} s [z,y], (3.11) 

which are visualized in Fig. lb). The sum on the right hand side of eq. (3.9) does not include nearest neigbours 
< x, y >, which are occupied by dimers or monomer pairs. Finally there exists the analogue of eq.(2.8) for the 
expansion of the spin-spin correlators projected on the spin S sector: 

l s (x, y)l s = S ^ x > y\ e K s)P{Ks). (3.12) 

AaS K S 
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IV. EXPECTATION VALUES OF THE BOLTZMANN FACTOR BETWEEN MONOMER- VALENCE 

BOND STATES 



Let us introduce the following notation for the eigenvalues hs and eigenstates \hs) of Hs'- 

H s \h s ) = h s V\h s ). (4.1) 

Furthermore let us consider the expectation values of the Boltzmann factor exp(/3Hs) between monomer- valence bond 
states \Ks): 

f(0,K s ) = (Ks\exp((3Hs)\K s )^Y, ex ^f 3Vh s)( h s\P(Ks)\hs). (4.2) 

h s 

Then eq.(3.9) -evaluated between the eigenstates \hs) - yields a system of linear differential equations describing the 
'evolution' of the expectation values /(/?, Ks): 

-^f((3, K S ) = (dV + 3N 1 (K S ) - N 2 (K s ))f((3, K s ) 

+ Yl (f(/3,Q+(x,y)K s )-f(p,Q_(x,y)K s )), (4.3) 

<x,y> 

with the initial condition: 

f(0 = O,K s ) = l, for all K s . (4.4) 

The expectation values f((3, Ks) contain the whole information on the thermodynamical and magnetic properties of 
the system: 

tr (exp/3tf s ) = — y2f(0,K s ). (4.5) 

aS Ks 

Eq.(4.5) results from the 'completeness' relation (3.4) for the projection operators P{Ks). 
Similarly one gets for the spin-spin correlators (x,y), weighted with the Boltzmann factor: 

tr [(x,y)exp(3H s ] —Vifcv] G K s )f(f3,K s ) (4.6) 
Aas k s 

as can be derived from the representation (3.12). 
In the combined limit: 

F,S^oo, M = S/V, (3 fixed (4.7) 
the free energy per site at fixed inverse temperature [3 and magnetization M is given by: 

F(f3,M) = ^logtr (expf3H s ). (4.8) 

Note that in the absence of an external magnetic field the magnetization M (in the antifcrromagnetic system) is zero. 
This means that the free energy F(/3, M — y = 0) in the spin S = sector contains already the whole information on 
the thermodynamics in the absence of an external field. Of course the same holds for the free energies in any of the 
sectors with fixed total spin S and for consistency they have to be the same. This can be seen immediately from the 
evolution equations (4.3), which differ for the various values of S by the number of monomers in the configurations 
Ks- It is at least plausible that a finite number of monomers will not change the behavior of the f((3,Ks) in the 
thermodynamical limit. 
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V. THE MONOMER DIMER CLUSTER EXPANSION 



In this section, we want to investigate an approximation scheme for the evolution equation (4.3). Let us first ask, 
what are the most important characteristics of the monomer valence bond configuration K$ which determine the 
variation of /(/3, Kg) with Kg. Obviously, such characteristic variables are the densities of dimers and monomer pairs 
x\ = Ni/V, x 2 = N 2 /V which appear explicitly in the evolution equation (4.3). It is clear from eq.(4.3) that the 
/(/3, Kg) increase much stronger with f3 on configurations with high dimer- and low monomer pair density. Next 
the clustering of the monomers and dimers seems to be important as will be demonstrated below. Therefore, the 
expectation values f([3,Kg) are supposed to depend on the numbers Nx(Ks), N 2 (Ks), ■■■ of dimers, monomer pairs, 
neighbouring dimer pairs, etc. 

Let us discuss first the lowest order in this approximative scheme, where we average the evolution equation (4.3) over 
all monomer valence bond configurations with fixed numbers Ni(Kg), N 2 (Kg) of dimers and monomer pairs. The 
approximation would be exact if the right hand side could be expressed as well in terms of the averages 

mNl ,N 2 ) = — 1 £ mK s ( Nl ,N 2 )). (5.1) 

V Ll M * Ks(N!,N 2 ) 

This is not possible for all the terms in the sum on the right hand side of eq.(4.3). Here we approximate the f(j3, Kg) 

he approximate evolution equation can be broi 

-f{P,N u N 2 ) = (dV + 3N 1 -N 2 )f(/3,N 1 ,N 2 ) 



by their average (5.1). The approximate evolution equation can be brought in the form: 

d_ 

0/3" 

+y E c h,n{Ni,N 2 , V)f(J3, iVi + h,N 2 + j 2 ). (5.2) 
31 ,h 

The coefficients Cj ± j 2 (Ni, N 2 , V) are given by : 

c^(NuN 2 ,V)= ] E E E 

v L > z ' <x,y> K(N 1 ,N 2 ) &= + ,— 



1=1,2 

Note that the permutations Q±(x,y) (eqs. (2. 13, 14), (3. 10, 11)) -applied on the monomer valence bond configuration 
Ks(N\) - can change the number of dimers and monomer pairs: 

N^Q^x^Ks) -N t [K s ) = ji±(K s ) = -2,-1,0,1,2 (5.4) 

N 2 (Q±(x, y)K s ) - N 2 (K S ) = j 2 ±(K s ) = -(d - 1), (d - 1). (5.5) 

The coefficients (5.3) are investigated in appendix B. They can be expressed in terms of probabilities to find certain 
clusters on the monomer valence bond configurations K$(Ni, N 2 ). In the combined limit: 

V, S, Ni,N 2 -> 00, xi = Ni/V, x 2 = N 2 /V, /3, M = S/V fixed, (5.6) 

they only depend on the densities x\, x 2 of dimers and monomerpairs: 

Cjij2(Ni,N 2 , V) -> c jlh (x u x 2 ). (5.7) 

Note, that the magnetization M — S/V - i.e. the monomer densities - is held fixed. It is a conserved quantity in the 
evolution process. 

The coefficients ( |5.3| ) acquire a particularly simple form in the limit: 

M — > 0, x 2 — > 0, c nn (xi 1 x 2 ) -> c n (xi)8 n oi (5.8) 

where 

c- 2 (xi) = -wx(xi), (5.9) 

c-i(xi) = 3iui(a;i) - w 2 (xi) - 2(2d- l)x t , (5.10) 

co(xi) — — 3wi(xi) + w 2 {x\) + 4(2d — l)xi — d + x\, (5-11) 

d(xi) = Wi(xi) + d- Xi(2(2d - 1) + 1). (5.12) 

Here wi(x),w 2 (x) are the probabilities to find on the valence bond configuration K$: 
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• dimer pairs with arbitrary orientation: w\{x\) 

• dimer pairs with parallel orientation: ^2(^1)- 

These types of clusters are shown in Fig. 2a-b for the case d = 2. 



VI. THE EVOLUTION EQUATION IN THE THERMODYNAMIC AL LIMIT 



In this section we are going to study the approximate evolution equation (|5.2j) in t he combined limit (5.6). Then 



(5.2) reduces to a partial differential equation for if we assume that the averages (|5.1|) have the form: 

/(/?, Ni,N 2 ) = exp(V<t>(J3, 11,3:2)), (6.1) 

J?-(f>(l3, x u x 2 )= log R t 1 = 1,2, (6.2) 

OXi 

d 

— <f)(f3, xx, x 2 ) = L(xi,Rx,X2, R2), (6.3) 



where 



The initial condition: 



L(xi, Ri,x 2 , R2) = d + 3xi - x 2 + ^ c Ji j 2 (x 1 ,x 2 )R{ 1 R? 2 ■ (6.4) 



Ox ,32 



(f>(/3 = 0, Xl ,x 2 )=0 (6.5) 
is a consequence of eq.(4.4). To solve eqs.(6.2-5) we first introduce -via a Legendre transform- a function: 

S(a:i, R u x 2 , R 2 ) = (j>[fi, x x , x 2 ) - p d< l>lP>*i> x *) (6 . 6) 

op 

which depends on the four variables x\ , x 2 , Ri , R 2 . The partial derivatives of $ with respect to these variables are 
not independent but obey the following four relations: 

"-(£-«**)(£)" ,=1 ' 2 (6 - 8) 

A solution of these equations is found with the ansatz: 

Ri R 2 

f f dR' dR' 

$(x 1 ,R 1 ,x 2 ,R 2 )=x 1 logR 1 + x 2 logR 2 - / / 1 2 g(R[,R' 2 ,L) 



R X R 2 

J^M,L)-J^ 92 (R 2 ,L), 



1 1 

Ri R 2 



(6.9) 



if the integrands g(R[, R' 2 , L), g±(R[, L), g 2 {R' 2 , L) with L = L(xi,R±, x 2 ,R 2 ) satisfy the conditions: 

R 2 / 

x l =x 1 (R l ,R 2 ,L)= I ^g{R 1 ,R' 2 ,L)+g 1 {R x ,L), (6.10) 
1 2 
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x 2 =x 2 (R u R 2 ,L) = J ^g(R[,R 2 ,L)+g 2 (R 2 ,L). 



(6.11) 



Note that the ansatz (6.9) satisfies the initial condition (6.5). According to (6.2), the initial point (3 = is mapped 
on the point R\ = R 2 = 1 in the Ri, R 2 -pl&nc. There the Legendre transform (6.6) has to vanish for all x\,x 2 . 
Eqs.(6.10,ll) define x\,x 2 as fuctions of R\,R 2 ,L with partial derivatives: 



Ri 



dxi 



dR 2 



Ri 



dx 2 



Ri,L 



= g(R u R 2 ,L). 



Therefore xi, I = 1,2 can be written as the gradient: 



x^—^—GiRuR^L) 
o log Ri 



(6.12) 



(6.13) 



of a 'potential' G(Ri, R 2 , L). Inserting (6.13) into (6.4) we end up with a partial differential equation for G(R\, R 2 , L). 
The solution of this equation demands the knowledge of the coefficients Cj 1 J2 {x\, x 2 ) in (6.4). They will be computed 
in the subsequent papers, where we treat the thermodynamical and magnetic properties of the antifcrromagnetic 
Heisenberg model in dimensions d— 1,2, 3. 



VII. SUMMARY AND PERSPECTIVES 



The main result of this paper is the evolution equation (4.3) for the expectation values (4.2) of the Boltzmann 
factor between monomer valence bond states. Once the evolution equation is solved, thermodynamical and magnetic 
properties (cf.(4.5)) as well as spin-spin correlators (cf.(4.6)) can be easily computed. The high temperature expan- 
sion of eq.(4.3) yields a recursion relation in a vector space of dimension vs{V). This is just the number of monomer 
valence bond configurations. This number can be reduced by using the symmetries of the Hamiltonian Hs- It seems 
to be realistic that a direct solution of such a recursion relation is feasible on clusters up to 18 sites. 
On the other hand, the exact evolution equation (4.3) can be considered as the starting point of a new type of approx- 
imation for the quantum Heisenberg model. In the thermodynamical limit the problem of solving the approximated 
cqation (5.2) can be reduced to the solution of a partial differential equation. The approximation of the evolution 
equation (4.3) can be systematically improved. The prize one has to pay is that the corresponding partial differential 
equations become more complex. 

We meet the most simple situation in the spin sector, where the approximated evolution equation (5.2) can be solved 
analytically- up to integrations. We will discuss this solution in a second paper and derive from it explicit expressions 
for the specific heat of the antifcrromagnetic Heisenberg model in dimensions d = 1,2,3. 
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APPENDIX A: THE PROPERTIES OF THE VALENCE BOND PROJECTION OPERATORS (VBPO) 



The properties of the VBPO's listed in eqs.(2.6), (2.8), (2.10) and (2.12) are based on the following anticommutation 
relations for the projection operators (x,y) defined in eq.(2.1): 

(x,yf = (x,y) (Al) 
l[(l,2),(2,3)] + = i((l ) 2) + (2,3)-(l,3)) (A2) 

i[(l, 2)(3, 4), (2, 3)]+ = 2)(3, 4) + (1, 4)(2, 3) - (1, 3)(2, 4)). (A3) 
For the proof of the completeness relation (2.6) it is sufficient to show that the sum over all VBPO's: 

^P(K)=^aj(Sy (A4) 

K J 

can be written as a power series in the square of the total spin S. Eq.(A4) obviously implies eq.(2.6): The right hand 
side -applied to any spin state- produces the same state. The left hand side -applied to any state with spin s > 0- 
yields zero. To verify (A4) we start with the sum: 

s i = ( Xl > yiX^, y2)-{x j ,y j ) (A5) 

xi,...,yj 

over the product of j disconnected valence bonds. We are going to show by complete induction that: 

Si = EM^ (A6) 

has the property we are looking for. Eq.(A6) is obviously correct for j = 1: 

Si = J2(x,y) =a w + ail S 2 . (A7) 

x=£y 

According to the definition (A5), the anticommutator [Ej,Ei]+ contains for types of contributions: 

• disconnected: The sites x, y do not overlap with the sites x\....yj. The sum over these contributions yields 

• one site connected: One of the sites x or y overlaps with the sites x\...yj. The sum over these contributions 
leads back to terms of the type £j, as can be verified with the anticommutator (A2). 

• two site connected: Both sites x and y overlap with xi,...,yj. Again the sum over these contributions leads 
back to terms of the type £j as can be seen from the anticommutator (A3). 

This completes the proof of eq.(A4). 

Next we turn to the proof of the representation (2.8). From eqs.(Al) and (A3) we get anticommutation relations for 
the VBPO's on two types of configurations: 

• K\. sites x and y are occupied by the same valence bond 

l[P(K 1 )Ax,y)}+ = P(K 1 ) (A8) 

• K 2 - sites x and y are occupied by two different valence bonds 

l -[P{K 2 ),{x,y)] + = ^(P(K 2 ) + P(Q + (x,y)K 2 )-P(Q_(x,y)K 2 )). (A9) 
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The action of the permutations Q±<j{x,y) 1 is defined in eqs.(2.13,14). From eq.(A9) we get 



[P(K 2 ) + P(Q_(x,y)K 2 ),(x,y)}+ = P(Q+(x,y)K 2 ). (A10) 
and by summation over all configurations of the type K 2 

J2lP(K), (x, y)} + (x, y)} + = ^ £ P{K X ). (All) 

K Ki Kx 

On the right hand side we have used the fact that the configuration Q + (x,y)K 2 is of the type K\. Evaluation of the 
left hand side by means of the completeness relation (2.6) and the projector property (A8) leads to the representation 
(2.8). Finally, the anticommutation relations (2.12) for the Hamiltonian Hq with the VBPO's are derived from 
cqs.(A8) and (A9) by summation over all nearest neighbour valence bonds (x,y). 

The properties (3.4), (3.9) and (3.12) of the projection operators P{K S ) can be derived in a similar way. For the 
proof of the completeness relation (3.4) we express the projection operators: 

s-i 

\[z u z 2S ] + ){[z u ...z 2S }+\ = c [] (,? 2 j(j + 1)) (A12) 

i=i 

as a power series of the square of the total spin s: 

2S 



s = 

1=1 



J2s(z t ) (A13) 



on the sublattice which is occupied by the monomers at sites zi, ...z 2 s- From the decomposition (A6) we know, that 
any power of s 2 : 
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(s 2 y = J2 c ^'i (Au) 

1=0 

can be written as a sum over terms EJ of the type (A5). In we sum over the product of I valence bonds on the 

sublattice with sites Zi, ...z 2 s- 

The sum over all projection operators P(Ks): 

iL v_ 

]T P{K S ) = ]T b^j = J2 ^0 2 y (A15) 

K s 3=0 j=0 

is then expressed in terms of the quantities Sj defined on the whole lattice. Eq.(A5) yields the expansion in powers 
of S 2 , which is sufficient for the proof of the completeness relation (3.4). 

For the proof of the anticommutation relation (3.9) and the representation (3.12) we have to generalize eqs.(A8) and 
(A9) to the corresponding equations for the projection operators P(Ks). Here, we have to consider three types of 
configurations: 

• K\. sites x and y are occupied by the same valence bond 

\[P{K 1 ),{x,y)] + =P{K 1 ). (A16) 

• K 2 : sites x and y are occupied either by two different valence bonds or by one valence bond and one monomer 

±[P(K 2 ),(x,y)] + = ±(P(K 2 ) + P(Q + (x,y)K 2 )-P(Q_(x,y)K 2 ). (A17) 

The action of the permutations Q±(x,y) is defined for the first case in (2.13,14) and for the second one in 
(3.10,11). 
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• K^: sites x and y are occupied by two monomers 

±[P(K 3 ),(x,y)} + = 0. (A18) 

From eq.(A17) we get 

[P(K 2 ) + P(Q-(x,y)K 2 ),(x,y)} + = P(Q+(x,y)K 2 ), (A19) 
and by summation over all configurations K 2 

J2l p (Ks), (x, y)}+- J2l p (Ki), (x, y)]+ ]T[P(X 3 ), (x, y)]+ (A20) 

K s Ki K 3 

Ki 

Finally, the representation (3.12) and the anticommutation relation (3.9) follow from (A16), (A18), (A20) and the 
completeness relation (3.4). 



APPENDIX B: MONOMER DIMER CLUSTERS ON THE CONFIGURATIONS K s 



We are going now to exploit equation (5.3) for the coefficients in the approximated evolution equation (5.2). The 
first sum on the right hand side of eq.(5.3) extends over those nearest neighbours < x, y > which are not occupied by 
a dimer or a monomer pair. The following situations - which are depicted for d = 2 and < x, y >=< 2,3 > in Fig. 
2a-f)-are left: 

• K a : dimer-dimer ( not parallel) 

• Ki,: dimer-dimer (parallel) 

• K c : dimer- valence bond 

• Kd- valence bond - valence bond 

• K e : dimer- monomer 

• Kf. valence bond- monomer 

The valence bonds in the situations K c , Kd, Kf must not be dimers in order to avoid double counting. On the 
configurations K a , Kb, K c , Kd the neighbouring sites < x,y >=< 2,3 > are occupied only by dimers and valence 
bonds. Therefore the permutations Q±(2,3) will not move the monomers on the configurations, which means that 
the number of monomer pairs does not change: 

N 2 (Q±(2,3)Kt) - N 2 (Kt) = for l = a,b,c,d. (Bl) 

The changes of the number of dimers (eq.(5.4)) on these configurations can be read of table I. 



TABLE I. 





K a 


K b 


K c 


Kd 


tf.(/*(l),A*(2),/i(3)) 


Kfi^x), ^(2)^(3)) 


31+ 


-1 








1 








31- 


-2 


-2 


-1 





-1 





32+ 











-1 


M(l)-M(3) 


fj,(x) - ^(3) 


32- 











-1 


M (2)-m(3)-1 


/x(2)- M (3)-1 
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Let us now turn to the configurations K e where monomers are involved. The permutations Q+(2, 3) and Q-(2, 3) 
will move the monomers from site 3 to site 1 and site 2, respectively - according to the rules given in Fig. lb. The 
resulting change (5.5) in the number of monomer pairs depends on the monomers surrounding the sites 1, 2, and 3 
on the configuration K e = K e {^{l), /z(2), /x(3)). If we denote by n(x),x = 1,2,3 the number of monomers on the 
neighbouring sites of x we find: 

N 2 (Q+{2, 3)K e ) - N 2 (K e ) = M (l) - M (3), 

7V 2 (Q_(2, 3)K e ) - N 2 {K e ) = M (2) - M (3) - 1. (B2) 

In the same way, one can see that the numbers [i(x), /i(2), /i(3) on the neighbouring sites of a;, 2,3 determine the 
change in the number of monomer pairs on the configuration Kf = Kf(/j,(x), /i(2), /ti(3)). 

The changes (5.4,5.5) in the numbers of dimers and monomer pairs - due to the action of the permutations Q+(2, 3) 
and Q_(2, 3) on the configuration Ki I = a, / - is summarized in table I. 

Next we need the probabilities wi (I = a, ...,/) to find these situations Ki on the monomer valence bond config- 
urations Ks(Ni 7 N 2 ). In principal, they can be determined by a Monte Carlo simulation of the monomer valence 
bond system. The situation simplifies considerably in the spin sector, where we have no monomers and where the 
complicated situations K e and Kf do not occur. The remaining situations K a — K4 have probabilities: 

w a = w 1 -w 2 , (B3) 
w b = w 2 , (B4) 
w c = 2(2d- l)ari - 2w u (B5) 
Wd = (d-xi)-w c -wi, (B6) 

(B7) 

which can be completely expressed in terms of the dimer pair probabilities wi(xi),w 2 (x 2 ) defined in section V. The 
coefficients (5.9-12) then follow from table I: 

C-2 = -w a - w b , (B8) 
c_i =w a - w c , (B9) 
co = Wb + w c - w d , (BIO) 
ci = w d . (Bll) 

(B12) 



FIGURE CAPTIONS 

1. The action of the permutations Q±(x, y) on a monomer-valence bond configuration K$'- 

• a) both sites x and y are occupied by valence bonds 

• b) site z is occupied by a monomer, site a; by a valence bond 

2. Different types of monomer-valence bond clusters at the nearest neighbour points 2 and 3. 



12 



yi 



x y 



x 2 



Q+(x,y)K s : 



Vi 



x y 



x 2 



X2 



Q-(x,y)K s : 



Fig. la) 



K s : z x y 
• 

z x y 
Q+(x,z)K s : 

z X 

Q-(x,z)K s : 



Fig. lb) 



K n : 



4 
3 



1 

2 

1 

2 



• 3 

4 

3 4 



K d : 



K,: o 



o o o 



1 2 3 



o o o 



Kf : o 



2 3 

Fig. 2a-f) 



